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SUMMAPy 


A low-power (1 S nN) laser velocimeter operating in the forward-scatter mode 
has been used to measure free-stream mean velocities in the Langley 0.3-Meter 
Transonic Cryogenic Tunnel. Velocity ranging from 51 to 235 m/s was measured 
with at least t1 -percent accuracy. These measurements were obtained for a 
variety of nominal tunnel conditions: Mach nun^rs from 0.20 to 0.77f total 

temperatures from 100 to 250 R* and pressures 101 to 152 kPa (1.0 to 1,5 atm). 
Particles were not injected to augment the existing Mie scattering material. 

It is postulated that the existing light scattering material in these tests was 
liquid nitrogen droplets normally injected to control the tunnel temperature. 
Signal levels obtained during the tests indicated that the average particulate 
diameter was greater than 1.0 m. Tunnel vibrations and thermal effects* which 
were considered to be potential problems before the tests* had no detrimental 
effects on the optical system. 


INTRODUCTION 

During the past several years* laser velocimeter (LV) techniques have been 
successfully applied to flow velocity measurements in aerodynamic test facili- 
ties (refs. 1 to 5) . The nonintrusive* LV technique reduces the probability 
of flow disturbances* often experienced with probe-type t.ichniques. In this 
technique* the velocity of particles present in the flow is equated to the gas 
velocity; thus the need is eliminated for calibrations or knowledge of other 
flow parameters such as temperature and pressure. 

With the increasing interest in aircraft efficiency at high speed* the 
aerodynamic flow characteristics for full-scale Reynolds number simulation in 
the transonic regime are import<mt. In a given wind tunnel with a specified 
gas* the Reynolds number can be increased by increasing the pressure and/or 
lowering the temperature. An accepted method for increasing the Reynolds 
number is by cryogenic cooling. 

This paper describes the results of LV tests conducted in the Langley 
0.3-Meter Transonic Cryogenic Tunnel. The purpose of the tests was to investi- 
gate the problems and potential use of this technique for measuring the flow 
velocity characteristics in this unique wind tunnel. A low-power one-component 
LV system operated in the forward-scatter mode was used to make a series of 
free-stream velocity measurements without a model in the test section. The 
naturally occurring Mie scattering particulates used in these tests are postu- 
lated to be residual droplets resulting from the injected liquid nitrogen (LN 2 ) 
which is used to maintain the desired temperature in the tunnel. The optics* 
electronics* data processing equipment and procedures* and general tunnel oper- 
ating procedure and characteristics are described in this paper. 

Also presented are results oC the LV tests which recommend that the laser 
velocimeter be further investigated for future use in t.iis facility. 



SYMBOLS 


d half-amplitude vibration displacement# cm 

E{ } expected value 

longitudinal function 
Fy vertical function 

f vibration frequency# Hz 

fg reference clock frequency# MHz 

fQ Doppler frequency# mz 

g gravitational acceleration 

M number of threshold crossings 

M free- stream Mach number 

OD 

N number of velocity measurements in ensemble 

n exponent in data from high-speed burst counter 

q = fo/Mfc 

Rjjy cross correlation function 

Sf fringe spacing# um 

S^y cross spectral density 

T integration time 

total temperature# K 

V absolute velocity component normal to fringes# nv^s 

Vjj calculated velocity, m/s 

individual velocity measurement# m/s 
Vm measured velocity# m/s 

X longitudinal direction 

Y true turbulence intensity, o/V 

Ym measured turbulence intensity 


2 



A uncertainty 

Cy error in mean velocity 

Cy error in turbulence intensity 

Hi randan variable 

0 angle between intersecting laser beams, deg 

X laser wavelength, run 

o true standard deviation, m/s 

0 |„ measured standard deviation, m/s 

T delay time 

A bar over a symbol denotes arithmetic mean. 


TECHNIQUE AND APPARATUS 
Wind Tunnel 


A schematic diagram of the Langley 0.3-Meter Transonic Cryogenic Tunnel is 
shown in figure ’ . It is a continuous-flow, fan-driven pressure tunnel which 
uses nitrogen as a test gas. The very cold temperatures required for high 
Reynolds number testing are achieved by direct injection of liquid nitrogen (LN2) 
just downstream of the test section. Stagnation temperatures from 77 to 339 K, 
stagnation pressures of 1.2 to 6.0 atm (1 atm « 1 01 kPa) , and Mach numbers from 
0.05 to 0.95 (the exact range depending upon the total temperature) can be 
attained by varying the rates of liquid nitrogen injection and gaseous nitrogen 
venting and the fan speed. The combined low temperature and high pressure can 
provide Reynolds numbers exceeding 330 million per meter. The tunnel is con- 
structed of tirelded 6061 -T6 aluminum alloy and encased in a fiberglass thermal 
blanket except for the test section access flange bolts and metal optics enclo- 
sures. The 12-bladed fan is driven by a 2.2-MW variable frequency motor. 


This facility is equipped with a 0.20 x 0.6I m, two-dimensional test 
section for airfoil testing which has four quartz windows. Two of these windows 
are located in the tunnel sidewalls and are D-shaped to permit observation of 
the flow over the top of an airfoil. The other two are pressure windows located 
on the seune optical axis between the plenum and the optics enclosures mounted 
on the test section. These pressure windows serve to contain the cryogen at 
whatever pressure the tunnel is being operated. 


The liquid nitrogen supply system consists of two, 106-m^ tanks which 
supply the liquid nitrogen to a pump. The pump supplies liquid nitrogen to the 
tunnel pressures up to 1.2 MPa (175 psig) at 0.0094 mVs to four digitally con- 
trolled LN2 supply valves. The liquid nitrogen is injected through these valves 


3 



downstream of the tunnel test section through four flush-mounted fan spray 
nozzles in the diffuser. The liquid nitrogen jet is atomized as it enters the 
tunnel flow. 

A venting subsystem allows removal of excess gas in order to maintain 
existing tunnel pressure as liquid nitrogen is injected to obtain the desired 
operating temperature. A more detailed description of tunnel design and opera- 
tion is given in reference 6. 


Laser Velocimeter Technique 

The basic principle of laser velocimetry depends on Mie scattering and the 
Doppler effect. The fundamentals involved can be explained by considering the 
fringe model shown in figure 2. The laser beam is split into two parallel beams 
of approximately equal intensity. A lens focuses the beams to a common point. 
Being coherent (temporal and spatial), the interacting wave fronts of the becun 
form interference fringes within an ellipsoidal region (sample volume) which 
is approximately centered at the focal plane. (The size of the sample region 
depends on the beam separation, the beam dieuneter, and the lens focal length.) 
The fringe pattern is a combination of alternating light and dark r«.gions of 
equal spacing. The spacing Sf is given by 


Sf 


X 

2 sin (6/2) 


where X is the wavelength of the laser light and 6 is the angle between 
the optical axes of the intersecting beams. If a particle in the moving fluid 
passes through the fringes, it scatters light from the two incident laser beams. 
This scattered light has a sinusoidal intensity variation with time. The fre- 
quency of intensity variation is a function of the fringe spacing and of the 
particle velocity in the direction normal to the fringes. The optical signal 
(Doppler shifted frequency) is converted by a photodetector into an electrical 
burst signal with frequency f^ given by 




2V sin (6/2) 


where V is the component of the absolute particle velocity perpendicular to 
the fringe pattern. 


Optical System 

The LV optics used in the tunnel tests discussed herein were configured 
to operate in the for ward- scat ter mode and to measure the axial component of 
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the flow velocity. A diagram of the system is shown in figure 3. It uses a 
He-Ne laser operating in the TEMqo n\ode and providing 1 5 nW of continuous power 
at a wavelength of 632.8 nm. The laser beam is divided into two parts by a 
50:50 beam splitter prism which yields a beam separation of 50 mm. A trans- 
mitter lens with focal length of 47.62 cm was used to focus the beams to a spot 
size (sample volume) 0.20 mm in diameter by 3 mm long. The light scattered by 
particles in the sample volume is collected by an FI :2 lens and is imaged on 
a high-gain^ low-noise photomultiplier tube (PMT) with an S-20 wavelength 
response by a 1 Ox (0.25 numerical aperture) microscope objective lens. The LV 
optical system was mounted in metal enclosures designed to house a schlieren 
system. This was done so that the enclosures could be purged with dry room 
temperature nitrogen and wrapped with heater tape to prevent moisture from 
condensing on the viewing windows (fig. 4). The optics installed in these 
enclosures are shown in figure 5. 


Signal Conditioning Electronics 

The signal conditioning electronics incorporated a 1 0-dB-gain wideband 
amplifier# a high-speed burst counter# and a minicomputer plus peripheral 
equipment. The high-speed burst counter is a real-time (time domain) device 
which uses a direct counting approach for measuring the time it takes for a 
particle to traverse a preset number of successive fringes. This time is 
measured with a resolution of 2 ns. The counter is equipped with band-pass 
filters for removing the pedestal voltage and with a double- threshold zero- 
crossing detector and comparator circuit for noise rejection. Data that have 
been validated are placed in a storage register which has a 16-bit binary out- 
put interface with a minicomputer. The data output format (10 data bits# four 
exponent bits) permits measuring a Doppler frequency range from 1.2 kHz to 
100 MHz. These data are fed into a minicomputer (binary form) which is pro- 
grammed to calculate the Doppler frequency from the following conversion 
equation: 


Frequency 


2.2 X lO’O 


2'' X Reading 


where the values of the exponent n and the reading are provided by the high- 
speed burst counter. For more detailed discussion of the digital counter, see 
reference 3. 


Data Acquisition and Processing 

The digital output from the high-speed burst counter is input to a 16-bit 
data acquisition and processing minicomputer where it is stored on a 0.5-megabyte 
floppy disk for later final processing. On-line data processing consists of the 
calculation of the basic statistical characteristics of an ensemble of velocity 
measurements. These results along with a plot of the measured velocity ensemble 
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histogram are displayed on a cathode-ray tube (CRT) terminal. The block diagram 
of the system is shown in figure 6. 

Data acquisition from the high-speed burst counter is controlled by the 
minicomputer through full handshake interaction. The maximum data acquisition 
rate of 20 000 measurements per second was sufficient to obtain the data in the 
present test since the average data rate was between 10 and 20 measurements per 
second. Each velocity ensemble consisted of 1024 measurements and the average 
data acquisition time was approximately 90 seconds. 

Once the velocity ensemble has been acquired, the ensemble arithmetic mean 
and standard deviation are calculated. Data measured by the high-speed burst 
counter greater than three standard deviations from the mean are considered 
spurious and therefore are deleted from the ensemble. The calculated mean value 
is then subtracted from the remaining data in the ensemble to maximize precision 
in the remaining statistical calculations, since single precision (6 decimal 
places) floating point numbers are used in the minicomputer. A new arithmetic 
mean and standard deviation are calculated based on the adjusted ensemble with 
the mean removed. If there are no spurious data points, the new mean is equal 
to zero; if there are, the new mean represents the difference between the pre- 
viously calculated mean and the true measurement mean. Thus the algebraic addi- 
tion of the previous arithmetic mean and the new arithmetic mean yield the true 
measurement mean velocity of the velocity ensemble. Likewise the new standard 
deviation represents the true standard deviation of the velocity ensemble since 
the spurious data were deleted from the ensemble. In addition to the mea'n and 
standard deviation, the following statistical quantities are also calculated: 
statistical uncertainty in the mean, statistical uncertainty in the standard 
deviation, normalized standard deviation, and the skew and excess of the 
ensemble. These results along with the histogrjup. of the ensemble (1.0 m/s 
resolution) were displayed on the CRT terminal and a hard copy of the CRT dis- 
play was produced. An example of the CRT display is shown in figure 7. 


ERROR ANALYSIS 

The accuracy of velocity measurement using the laser velocimeter is depen- 
dent upon the accuracy with which the particle follows the fluid flow, gradients 
within the flow field, errors derived from the optical and electronic systems, 
and measurement uncertainties due to the statistics of the measured velocity 
ensemble. In the present test, the free-stream velocity was measured in an 
empty test section; there were no detectable effects on the measurements due 
to velocity gradients within tlie velocimeter seunplo volume. Further, the 
velocity lag experienced by the particles as the flow accelerates in the tunnel 
throat was ,yot apparently significant since the standard deviations of the 
measurement ensembles were small and reasonably symmetric, even though the dis- 
tribution oL the liquid nitrogen particle sizes was suspected to be large 
(»1.0 Urn in diameter) based on the signal levels received from the laser 
velocimeter. Since the velocity gradients were negligible and the measured 
standard deviations were small, the following errors should also be negligible: 
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1 . Particle tracking errors (inability of particles to follow the flow 
field and related velocity broadening due to the size distribution of the 
particles) 

2. Velocity bias (the probability that a larger number of faster particles 
will be measured than slower particles per unit time) 

3. Measurement dead zone (particles passing through the sample volume with 
a trajectory that crosses too few fringes to yield a measurement) 

4. Velocity gradient errors (velocity gradient across the sample volume 
yielding a broadening of the velocity distribution which is not due to turbu- 
lence) 

5. Sample volume position (unfaithful representation of the flow field due 
to position inaccuracy) 

However, the remaining error sources must be considered; 

1 . Cross-beam angle measurement error and nonpar all.l ism of the fringes 

2. Threshold limit error, clock synchronization error, and quantizing 
error in the high-speed burst counter 

3. Uncertainties in the statistical calculations of the measurement 
ensembles 

The two error sources due to the optical system are cross-beam angle 
measurement and fringe nonparallelism. The cross-beam angle error is a velocity 
bias due ta the uncertainty in locating the center of the laser beam during the 
geometric measurement of the angle. The uncertainty band is estimated to yield 
a bias uncertainty in velocity of tO.6 percent for the present optical system. 

The fringes within the sample volume diverge if the spatial location of the waist 
of each laser beam is not located at the beam crossover point. In the present 
system, the waist for each beam extends 3 cm on both sides of the crossover 
point; thus, the errors due to fringe nonparallelism are reduced to an insig- 
nificant level. 

The three error sources in the electronic systems, threshold limit error, 
clock synchronization error, and quantizing error, are developed in detail in 
the appendix. Although the threshold limit error does not apply for the high- 
speed burst counter used in the present test, the clock synchronization and 
quantizing errors do and are dependent upon the mean velocity. As an example 
of their magnitudes, a worst case condition is used: mean velocity of 235 ia/s 

and an ensemble standard deviation of 1.68 m/s. Clock synchronization yields a 
bias error in the mean velocity of 0.47 percent and a bias error in the 
normalized standard deviation of 0.47 percent. Quantizing yields a random 
error in the mean velocity of 0.0003 percent and a bias error in the normalized 
standard deviation of 29.4 percent. In the present test, the normalized stan- 
dard deviation measurements were corrected for quantizing error, but not for 
clock synchronization error. The mean velocity bias error due to clock 
synchronization was not removed. 
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Estimates of the system eccoc sources in the velocity are summarized 
in the following table; 


Velocity error source 

Bias, 

percent 

Random, 

percent 

Cross-beam angle 

±0.6 


Diverging fringes 

Negligible («0.01) 


Clock synchronization 

0.47 


Quantizing 


Negligible (<<0.001 ) 


The final error source is the statistical uncertainty in estimating the 
true mean and standard deviation of the velocity of the flow. These errors 
depend on the number of particle velocity measurements included in the ensemble 
and the magnitude of the standard deviation of the ensemble. These uncertain- 
ties, presented in table I for each ensemble, average approximately 0.03 percent 
for the mean velocity and 3.5 percent for the standard deviation. 


TEST PROCEDURE 

The tunnel is cooled prior to operation by running the drive motor and fan 
at a low speed and injecting liquid nitrogen at a low flow rate. After reach- 
ing the desired operating temperature, the tunnel flow is brought to the 
desired Mach iiumber and pressure by adjusting the drive fan speed, controlling 
the quantity of liquid nitrogen being injected, and regulating the amount of 
nitrogen gas vented to the atmosphere. For the tests discussed herein, the 
nominal free-stre{im temperatures range frcm 100 to 250 K and Mach numbers frcm 
0.20 to 0.77. The laser velocimeter was adjusted to measure the axial free- 
stream velocity component along the test section centerline as shown in figure 8. 
The test section was empty. For these tests, naturally occurring particulates 
(postulated to be LN2 droplets) provided the light scattering particles neces- 
sary for the LV operation. Vibration levels on both sides of the test section 
were measured with a triaxial accelerometer mounted on the optics enclosures 
approximately 0.5 m from the test section wall. 


RESULTS AND DISCUSSION 

While preparing the tunnel for this series of tests, edge cracks were 
found in one of the pressure windows. Safety considerations required that the 
tunnel be operated at the lowest practical differential pressure without exceed- 
ing 1.0 to 1.5 atm in the plenum. Therefore, the conditions were limited to 
a Mach number range from 0.20 to 0.77 and a stagnation pressure of 1.5 atm. 
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At the colder operating temperatures (TOO K), the tunnel test section walls, 
with the viewing ports, could contract by approximately 0.32 cm. Refraction 
effects caused by this dimensional change had no detectable effect on the bV 
optical alignment. Thermostatically controlled heater tape wrapped around the 
optics enclosures kept the internal ambient temperature above 266 K. Dry nitro- 
gen gas purged through the optics enclosures kept the viewing ports free of 
condensation. 


Vibration Measurements 

The induced vibrations most critical to the optical alignment are the 
longitudinal mode (parallel to the flow) and the vertical mode (normal to the 
flow), with the worst cases occurring at the higher Mach numbers and tempera- 
tures. Typical results of vibration measurements for the transmitter and 
receiver optics enclosures (figs. 9 and 10) show distinct peaks at 100 and 400 
Hz. Since the displacement varies ir versely with the square of the frequency, 
100 Hz represents the larger displacement. At 100 Hz, the longitudinal and 
vertical vibration levels at the transmitter and receiver were on the order 
of 0.32g. The half-amplitude displacement, approximated by 


d 


24.84g 

f2 


is 8 X 10“^ cm* Coupling between the vibration measurements for the transmitter 
and receiver enclosures was determined by computing the cross spectral densityr 
which is expressed as 


Rxy(t) di 


Sxy 


I 


where lU is the angular frequency, j ■ |T-1 , and Rxyd) is the cross corre- 
lation function given by 


lim 



Fx(t) 


Fy(t - T) dt 


with T being the integration time and t the delay. Typical cross spectral 
density plots are shown in figure 11 . For the range of test conditions, the 
vibration levels were not severe enough to affect the optical aligrbsent or the 
laser operation. 
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Velocity Measurements 


Free-stceam axial mean velocity measurements as a function of Mach number 
and temperature are Lhown in figure 12. A comparison of the measured '''locities 
;v,n) ®nd those calculated from tunnel parameters (V^.) (sho»*n in tig. 1 J) indi- 
cates good agreement, typically better than 1 percent. The standard deviations 
for this plot cannot be shown because they would appear within the test-point 
symbols. A plot of the normalized standard deviation is shown in fig- 

ure 14. Over the test Mach number and temperature range, the normalized stan- 
dard deviation was usually less than 1 percent. The velocity measurement data 
and nominal tunnel conditions are listed in table I. The size of the particles 
(Il<2 droplets) from which these measurements were obtained was not determined, 
although it was possible to estimate frcm the burst patterns observed on an 
oscilloscope that many of the particles were large (»1.0 urn in diameter). 
Although particle size is usually important to laser velocimetry accuracy, it 
was not found to be critical in this instance because only the free-stream 
velocity was measured without a model in the test section and velocity gradients 
were small. The typical rate of occurrence and concentration of the measured 
particles were on the order of 10 to 20 per second and 10^ per cubic meter, 
respectively. 


CONCLUSION 

Free-stream velocity measurements have been made in the Langley 0. 3-Meter 
Transonic Cryogenic Tunnel using a low-power (15 nM) laser velocimeter. These 
measurements were obtained in an empty test section for the foMowing nominal 
tunnel conditions: Mach numbers from 0.20 to 0.77, temperatures from 100 

to 250 K, and pressures 1.0 to 1.5 atm. The laser velocimeter measurements 
typically agreed to within 1 percent with velocities calculated from tunnel 
parameters. The overall normalized standard deviation ((^n/Vni) 

1 percent. Vibrations and the inherent tunnel contraction due to temperature 
changes had negligible effect on the optical alignment and the laser operation. 
These favorable results recommend that laser velocimeter be further investigated 
for future use in this facility. 


Langley Research Center 

National Aeronautics and Space A<kiinistration 
Hampton, VA 23665 
March 10, 1981 
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APPENDIX 


CORRECTION FOR THE EFFECTS OF SIGNAL PROCESSING BY HIGH-SPEED 
BURST COUNTER ON TURBULENCE INTENSITY MEASUREMENTS 

by 

James F. Meyers 

The high-speed burst counter measures the time for a fixed number of cycles 
in the Doppler burst; errors occur in the normalized standard deviation (or 
turbulence intensity measurement) due to the digitizing of this time by a digi- 
tal reference clock. This effect is negligible until the flow turbulence inten- 
sity approaches 1 perc?n*c# sinc.e the true velocity probability density function 
is described by a histogram of only a few velocity increments. The limit of 
this effect is when the probability density function is described by a histogram 
containing only two increments, e.g., 625 and 626 reference clock counts. 

In order to determine the equation necessary to correct the data for this 
effect, the error equations must be developed. In general, there are three 
error sources due to the high-speed burst counter: (1 ) threshold limit error, 

(2) clock synchronization error, and (3) quantizing error. The threshold limit 
error is due to the use of a nonzero volt Schmidt trigger to digitize the 
Doppler burst. Since the high-pass filtered Doppler burst has a Gaussian 
envelope, a constant voltage trigger trips at different phases for adjacent 
cycles in the Doppler burst. Present state-of-the-art high-speed burst couplers 
use a true-zero-crossing detector which eliminates this error if the Doppler 
signal is properly filtered. Since tnis was the case in the present test, thi5"> 
error source is not discussed further. The clock synchronization error is due 
to the mismatch of the randomly occurring Doppler burst with V\e start of a 
clock cycle. Th« result of this error is a 1 -clock-count uncertainty in the 
measurement for Doppler bursts containing signals of the same frequency. The 
quantizing error is due to the digitizing of time during which a fixed number 
of Doppler cycles occur. Since the time to be measured is a continuous analog 
function, the quantizing error has an uncertainty of clock count. 

Once the aforementioned error sou«,ces have been described, they may be 
combined for a total error in the mean velocity measurement and a total error 
in the turbulence intensity measurement. From these equations, the error cor- 
rection equations may be developed and used to correct the measurements made 
with the laser velocimeter. 


Basic Statistical Equations 

The ^wo quantities in the error analysis described herein are the mean 
and the standard deviation of the measured velocity ensemble. The standard 
deviation measurement is then used to solve for the turbulence intensity y 
( equal to o/V) of the flow. The expected value of the measured velocity (or 
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mean velocity, Vnj) may be approximated by a Taylor series in terms of the 
true mean velocity V; 


e{v„} = Vm = f(C) + (X - C)f'(?) + 


'm 


2 ! 


f”(0 + 


(Al) 


where 


f {^) = E{Vi) = V 

= e{V| 2} = v2(1 + Y^) 
f " {^) = E{Vi^} = V3(l + 3Y^) 


and where the expected value of the random variable is 


(X - = ECrii^} 

Note that since and are independent, 

(X - = ECriiVi^} = E{ni) E{Vi2} 


Thus, 


V„, = V + v2(1 + Y^) Eirii) 


V^(l + 3Y^) 

+ 

2 


E{ni2} 


The definition of variance (o«2) is 


0 2 = 2} Y 2 


(A2) 


where 


£{7^2} = E{(Vi + t1iVi2 + rii2Vi3)2} (A3) 

and since and are independent, 
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= E{Vi2} + 2 E{ni> E{Vi3} + 3 £{11^2} E{Vi^} 

+ 2 E(ni^} E{Vi5) + E{rii^} E{Vi6} (A4) 

The moments of Vj are 

E{Vi2} = v2(1 + Y^) 

EWi^} = v2(l + 3 y^) 

E{Vi^} = V^(T + 6 y^ + 3y^) 

EfVi^} = v5(1 + IOy^ + 15y^) 

E{Vi®} =0(1 + 15 y^ + 45y^ + 5y®) 


T‘ is, 

E{V„2} = 0(T + y2) + 20(1 + 3y^) E{r«i} + 30(1 + + 3y4) Eihi^} 

+ 20(1 + IOy^ + ISy'*) E{rii^} + 0(1 + 1 5 y^ + 45y^ + 5y^) E{ni^} (A5) 

and from equation (A2 ) , 

V„2 = 0 + 20(1 + Y^) E{r)i} + 20(1 + 3 y 2) E{rii2} + 0(1 + y^) 2 E2{rii} 

+ 20(1 + Y^) (1 + 3y^) Elm} E{ni^} + 0(1 + 3 y^) 2 E2{rii2} (A6) 


Then, 

o„2 » o2 + 40y2 E{rii} +0(1 + 12 y^ + 9y^) £{^2} - 0(1 + y^)^ E2{rii} 

+ 20(1 + IOy^ + 15y^) E{ni^} - 2 v5(1 + y^) (1 + 3y 2) E{ni) E{ni^} 

+ 0(1 + ISY^ + 45y^ + 50) E{rii^} - 0(1 + 3 y 2)2 E2{rii2} (A?) 
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appendix 


Clock Synchronization 

The uncertainty in the clock synchronization is a uniform random variable 
within the band from 0 to -1 clock pulse. Thus the functions of become 


wher 


E{nil = V ^ 


1 


2\MSffci 2V\Mf^/ 2V 


E{rii2} 


V— f = ^ 

3\MSffc/ 3v2\Mfc/ 3 v2 


Elhi^} = V ’ 


4\MSffci ' 4Wfc; ' 4 v3 


E{rii4} 


5\MSffc/ ~ 5v4\Mfc/ ' 5 v4 


M 


number of Doppler cycles used 


Sf fringe spacing 

reference clock frequency 




q 


V 


Doppler frequency = — 

Sf 


Mfc 


In the case under study, the number of Doppler cycles used is eight, the 
reference clock frequency is 500 MHz, and the maximum Doppler frequency possible 
is 100 MHz. The maximum value of q is 0.025 which is much greater than q^; 
thus all terms q*^ with n S 2 in equations (A2) and (A7) are neglected. 
Equations (A2) and (A7) now become 
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V„ = V + v2(1 + y2) E{rii} = V 

= o2 + 4V^Y^ E{ni) = 0 ^ (1 + 2q) 


1 + -n + Y^) 
2 


The error in the mean velocity Cy is 


Vm - V 


£v = 


(1 + Y^) 


(A8) 

(A9) 


(AlO) 


and the error in turbulence intensity £y is 


= 


Vm Y q 

• -(1 - Y^) 

Y 2 


(All) 


Quantizing Error 

The uncertainty in the timing of the Doppler signal due to the digitizing 
of time by the reference clock is a uniform random variable within a band from 
-1 clock pulse to +1 clock pulse. Thus the functions of Hi become 


E{rii} = 0 


E{r)i2) 


1/ 1 \‘ 1 /fn r q2 


3\MSf fc/ 3v2\Mfc I 3v2 


E{ni2} = 0 




-f— V “ — 1 = — 

5\^MSffc/ 5v4\Mfc/ 5 v4 


Then equations (A2) and (A7) become (q^ terms neglected) 
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Vm 


V + v3(1 


+ 3y2) 


/qM _ 

q^ 2 

> 

II 

l> 

1 

1 + — (1 + 3Y^) 
3 


(A12) 


0^2 = o2 + v^(1 


+ 12y2 + 9Y^) 



o2 


v2q2 
+ (1 


3 


+ 1 2y^ + 9Y^) 


(A13) 


The error in the mean velocity Cy is 


Vm - V q2 

Cy = — r = — (1 + 3Y^) 

V 3 


(A14) 


and the error in turbulence intensity Gy 



(A15) 


Data Correction 

To correct the data on the basis of the equations derived for the afore- 
mentioned high-speed burst counter errors, the nature of the errors must be 
determined and their significance understood. The error in the mean velocity 
measur ,*inent due to clock synchronization is biased toward a higher velocity 
than the true velocity since the measurements are always correct for 1 clock 
count less. The error in the mean velocity measurement due to quantizing error 
is a random error since the probability that there wi'l be one too many clock 
pulses is equal to the probability that there will be one too few. Both errors 
in the turbulence intensity measurement are positive bla*^ errors. 

The total mean velocity error is determined by combining equations (Al 0) 
and (A14' 


= -(1 + ± — (1 + 3y2) (A16) 

2 3 


However, as stated above, q is much greater than q^; thus the second term 
is neg’ected and only a bias error described by equation (Al 0) is considered. 
Base'* on equation (A10), the true mean velocity nay be found from 


U 





APPENDIX 


2V„ 

V = 

2 + q(1 + y2) 


(AT 7) 


The turbulence intensity error is determined by combining equations (All) 
and (A15): 





- v2 


Y^) 



(A18) 


At large values of turbulence intensity^ the first term is dominant even though 
its magnitude is very small. However as the turbulence intensity decreases, 
the second term becomes dominant and its magnitude becomes significant. Thus 
the error correction is developed from the second term, which is the quantizing 
error given in equation (A1 5) , and used to correct the data at low turbulence 
intensities. The resulting expression is then simplified using binomial series 
expansion techniques which yield the following expression: 


Y - 



(A19) 


In order to verify equation (A19), an extremely stable FM generator was 
used as a single source for a high-speed burst counter and a frequency tracker. 
The output from the frequency tracker was used as the standard for determining 
the input turbulence intensity. The measurement errors for the high-speed 
burst counter as a function of input turbulence intensity are shown by the 
symbols in figure 15. The curve is the expected error based on the calculations 
from equation (Al 5) . Figure 16 contains the results from comparing the measured 
turbulence intensity from the high-speed burst counter with the input turbulence 
intensity (circular symbols). If the measured data are correct, they lie along 
the 45^ line. As shown in the figure, the quantizing error becomes significant 
at a turbulence intensity of 1 percent. When the turbulence intensity decreases 
to the point where all the data are contained within two reference clock counts 
(the horizontal line at 0.36 percent for this test case), the data and any cor- 
rection of the data arc unreliable. By applying equation (A19), the data are 
corrected and are shown by the triangular symbols. 

Thus the errors in the velocity measurements due to the high-speed burst 
counter may be removed since both the error in the mean velocity measurement 
and the error in the turbulence intensity measurement are bias errors and may 
be removed using equations (A17) and (A19). 


\ 
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Figure 3.- Diagram of the laser velocimeter system. 
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Figure 4 .- Metal enclosure for optical system mounted on 
tunnel test section. 


23 














FILE NUMBER 


♦ 


>1 

ns 

IH 

04 

tn 

n 


•*-4 

*0 

(U 

o 

(U 

ns 

4J 

r-i 

nj 


w 

03 


0) 

Tl 


r- 

♦ 

c 

nj 



§ 

u 


u> ^ 

■ 

(0 (ft 01 

D> 


(ft Vn (ft 

0 

jj 

(ft 

V © ^ 

U) 

f- 

O (ft • (ft CD 

X7 

z 

• *T Ift • 


M 

• ft" 

fe 

O 

> ♦ 

CZi 

U 

QL 

1- 



M 

r-f 

ns 

<C 

o 

0 

1- 

o • 

‘H 

CL 

<n 

•-la a 

V*4 

i: 

o 

UJ a <Z 


Z>Z E a (ft 

1 

u. 

<L<ZO (ft 

• 

o 

ZUlEM3y 

r- 


CEO(fty O 

a> 

D 

W M X 


Z 

E^(fto-(fty 

mJ 



•H 

fra 


PLENUM 


i 


P1X)R 


-> , 


;tv 



00 

0 

1 

0 
00 

1 



27 




AUTO POWER SPECTRAL OENSIT 
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(a) Transmitter enclosure. 
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(b) Receiver enclosure. 

Figure 9.- Acceleration power spectral density for induced 
vertical-mode vibrations of LV optics enclosures. 

M - 0.77; Tt « 250 K. 
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(b) Receiver enclosure. 

Figure 10.- Acceleration power spectral density for induced 
longitudinal'inode vibrations of LV optics enclosures. 

« 0.77; Tk « 250 K. 
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Figure 13.- Measured mean velocity versus calculated mean velocity 
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